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Abstract. Pseudo-differential and Fourier series operators on the torus T" = 
(]R/27rZ)" are analyzed by using global representations by Fourier series instead of 
' local representations in coordinate charts. Toroidal symbols are investigated and 

, the correspondence between toroidal and Euclidean symbols of pseudo-differential 

' operators is established. Periodization of operators and hyperbolic partial differen- 

tial equations is discussed. Fourier series operators, which are analogues of Fourier 
integral operators on the torus, are introduced, and formulae for their compositions 
, with pseudo-differential operators are derived. It is shown that pseudo-differential 

and Fourier series operators are bounded on under certain conditions on their 
, phases and amplitudes. 
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^ ■ 1. Introduction 

In this paper we investigate a global quantization of pseudo-differential operators on 
the torus. For this, we develop necessary elements of the toroidal microlocal analysis. 
Using the toroidal Fourier transform (Fourier series) we can quantize the usual classes 
of pseudo-differential operators yielding a full symbol of an operator on the torus. 
The main difference with the (Kohn-Nirenberg) quantization in Hormander's symbol 
classes ^^^(T" x M") (which we will often call Euclidean quantization to emphasize 
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the difference) is that the toroidal symbol belongs to the corresponding symbol class 
S^giT"^ X Z") with frequency variable ^ G in the integer lattice; this quantization 
will be called the toroidal quantization. We will analyze the relation between these 
symbol classes and between corresponding pseudo-differential operators. 

We will also investigate the corresponding toroidal version of Fourier integral oper- 
ators. To distinguish them from those defined using the Euclidean Fourier transform, 
we will call them Fourier series operators. The use of Fourier series will allow us to 
obtain a global representation of these operators, thus making trivial a number of 
topological obstructions known in the standard theory of Fourier integral operators 
on manifolds. We will prepare the machinery and describe how it can be further 
used in the calculus of Fourier series operators and applications to hyperbolic partial 
differential equations. In fact, the form of the required discrete calculus is not a-priori 
clear. For example, composition formulae for pseudo-differential operators rely on the 
Taylor expansion of symbols which are now defined on the discrete lattice Z". Thus, 
we develop the corresponding versions of the periodic and discrete analysis which are 
still quite similar to formulations in the standard Euclidean theory. In particular, this 
will include the analysis of differences and some elements of the microfocal analysis 
such as toroidal wave front sets, etc. We study pseudo-differential operators in detail 
giving an explicit relation between Euclidean and toroidal symbols. Moreover, we 
will describe how this relation can be used to relate Euclidean and toroidal quanti- 
zations of the the same operators using a periodization operator that we analyze for 
this purpose. 

However, our analysis will reveal certain limitations for the development of an 
equivalent full theory for Fourier series operators. The main difference with pseudo- 
differential operators is in the behaviour of the wave front set under the action of 
operators. In the case of an pseudo-differential operator, the wave front set does not 
move and so we remain in the space- frequency space T" x Z*^. In the case of a Fourier 
series operator, their integral kernel may have singularities away from the diagonal 
in T" X T*^, which means that under the flow the wave front set may no longer be 
a transformation of T" x Z^. Nevertheless, we will be able to resolve this issue in 
the composition formulae by using (unique modulo smoothing) extensions of toroidal 
symbols. 

It was realized already in the 1970s that on the torus one can study pseudo- 
differential operators globally using Fourier series expansions, in analogy to the Eu- 
clidean pseudo-differential calculus. These periodic pseudo-differential operators were 
treated e.g. by Agranovich [H E]. Contributions have been made by many authors, 
e.g. by Amosov [Sj, Elschner [TT], McLean [16], Melo [17], Prossdorf and Schneider 
|18j . Saranen and Wendland [26], Turunen [31], Turunen and Vainikko [32j, Vainikko 
and Lifanov [331 El], and others. However, most of these papers deal with one- 
dimensional cases or with classes of operators rather than with classes of symbols. 
In this paper we in particular develop the symbolic analysis of pseudo-differential 
operators on the torus thus closely relating Euclidean and toroidal quantizations of 
the same classes of operators. 

We note that the use of operators which are discrete in the frequency variable allows 
one to weaken regularity assumptions on symbols with respect to ^-variables. Symbols 
with low regularity in x have been under intensive study for many years, e.g. see 
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Kumano-go and Nagase [15j, Sugimoto [29J, Boulkhemair Garello and Morando 
|12j . and many other contributions. However, in these papers one assumes symbols 
to be smooth or sufficiently regular in ^. The discrete approach in this paper will 
allow us to reduce or to completely remove the regularity assumptions with respect 
to ^. For example, no regularity with respect to ^ is assumed for estimates, and 
for elements of the calculus. We note that a notion of pseudo-differential operators 
with low regularity in ^ can still be introduced, as was done e.g. by Sjostrand [2Z], or 
with symbols in more general modulation spaces, see e.g. Grochenig [T3j. Although 
such operators form an algebra, no explicit symbolic calculus of such operators is 
available. In this respect the toroidal quantization has an advantage of looking only 
at discrete values of ^, thus removing the regularity issue altogether, and yielding the 
symbolic calculus. 

We note that the approach developed in this paper may be generalised to other 
(non-commutative) Lie groups. Thus, we will develop global quantizations of pseudo- 
differential operators on in [23] by identifying it with group SU(2) via the quater- 
nionic product in M^. Moreover, in [21] we will also develop global quantizations of 
pseudo-differential operators on arbitrary compact Lie groups and on certain homoge- 
neous spaces without resorting to local coordinate charts. We note that this analysis 
is rather different from the one developed by Taylor in [30] because of the obtained 
descriptions of symbol classes and their transformations. In any case, the relation 
between compact Lie groups and M" is most transparent in the case of the torus and 
the utilisation of this relation brings considerable simplifications in the analysis of 
this paper compared to that in [23j and [24J , which rely heavily on the representation 
theory. Some results proved in this paper were announced in [22]- 

We fix the notation for the torus as T" = (]R/27rZ)". Often we may identify T" 
with the cube [0,27r)" C M"" (or [— 7r,7r)"), where we identify the measure on the 
torus with the restriction of the Euclidean measure on the cube. Functions on T" can 
be thought as those functions on MJ^ that are 27r-periodic in each of the coordinate 
directions. We will often simply say that such functions are 27r-periodic (instead of 
27rZ'^-periodic). 

For e N we will also use the notation A^T" = (M/27rA^Z)" and call this the 
N -inflated torus, or simply an inflated torus if the value of A^ is not of importance. 
One of the reasons to use large A^ in this analysis is that if one wants to embed a 
compactly supported problem in M" into a torus, one may need to choose A^ large 
enough for technical simplifications. We will briefly discuss these issues. 

One particular application discussed in this paper will be to hyperbolic equations on 
the torus. In particular, the developed analysis can be applied if we embed certain 
problems in M" into the torus. In general, we can observe that due to the finite 
propagation speed of singularities in hyperbolic equations we can usually cut-off the 
equation and the Cauchy data for large x in the local analysis of singularities of 
solutions for bounded times. Then we can embed the problem into T" using rescaling 
(or into the A^-inflated torus A^T" without rescaling) in order to apply the periodic 
analysis developed in this paper. One of the advantages of this procedure is that since 
phases and amplitudes now are evaluated only at ^ G Z"" one can apply this also for 
problems with low regularity in i^. This type of application is particularly important 
for weakly hyperbolic equations or systems with variable multiplicities when at points 
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of multiplicities characteristics become irregular. If the principal part has constant 
coefficients then the loss of regularity occurs only in ^, so techniques developed in 
this paper can be applied. Such applications to systems will appear elsewhere and in 
this paper we shall briefly discuss only the case of scalar equations. 

We also note that this approach to the representation of solutions is related to 
the work of Bourgain [HI El E] on Strichartz estimates for solutions to Schrodinger 
equations on the torus T". Bourgain used the Fourier series representation 

u{t,x)= ^'^fix)=J2 e'^^-«-*l«l')/(0, 

which shows in addition that the solution is periodic in time. This representation 
allowed him to prove, for example, the equality ||w||^4(']pxT") ~ II^^IIl2(txT") leading 
to the corresponding Strichartz estimates and global well-posedness results for non- 
linear equations (while dispersive estimates fail even locally in time). In Section [TO] 
we discuss solutions to hyperbolic equations with variable coefficients but leaving 
corresponding nonhnear applications outside the scope of this paper. 

Overall, the analysis developed in this paper lays down a foundation for tackling a 
variety of global problems on the torus with the use of the full toroidal symbols (ques- 
tions like global hypoellipticity, global solvability, estimates for pseudo-differential 
operators, Fefferman-Phong inequality, etc.), and for applying these techniques to 
partial differential equations. Again, we leave further applications outside the scope 
of this paper. 

In Section[2]we fix the notation used throughout the paper. In Section[3]we develop 
the necessary machinery for the discrete calculus, in particular we prove the Taylor 
expansion formula on the lattice giving an estimate for the remainder. Section 
m is devoted to the notion and main properties of toroidal symbols together with 
formulae for their calculus. In Section [3] we establish the relation between Euclidean 
and toroidal symbols and in Section [6] the relation between operators with Euclidean 
and toroidal quantizations. In Section [7] we prove basic results about toroidal wave 
front sets. Section M introduces Fourier series operators and estabhshes composition 
formulae between Fourier series and pseudo-differential operators. In Section [9] we 
prove results on the L^-boundedness of pseudo-differential and Fourier series opera- 
tors. Finally, Section [10] is devoted to an application to hyperbolic equations showing 
a way how to embed such equations into the torus by periodization. 

We will use the following notation in the paper. Triangles A and A will denote 
the forward and backward difference operators, respectively. The Laplacian will be 
denoted by =Sf. The Dirac delta at x will be denoted by 6x{y) and the Kronecker 
delta at ^ will be denoted by 5^,^. We will use the standard notation for multi-indices 
a e Nq, with No = {0} U N, for example denoting |a| := X]j=i — ' ' '(^n-: 

= d^l ■ ■ -5"", etc. We will also use the standard notation D^^ = —idx^ = y^f-, 

where i = a/— 1. We will write a G (7°°(T" x Z") when function a(-,^) is smooth on 
T" for all ^ G Z". For symbol classes S^^, we will often write simply S"^ for the class 
S]^Q. For ^ G M", we will write (^) = (1 -|- I^P)"*^^^. To avoid keeping track of constants 
in integrals we will normalise the measure by defining dx = (27r)~"' dx. Constants 
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will be usually denoted by C (sometimes with subscripts) , and their values may differ 
on different occasions, even when appearing in subsequent estimates. 

2. Preliminaries 

In this section we will introduce some notation which will be frequently used in the 
sequel. Let iS(M") denote the space of the Schwartz test functions, and let iS'(R"^) 
be its dual, i.e. the space of the tempered distributions. The Dime delta eomb 
Szn : 5(M") ^ C is defined by 

and the sum here is absolutely convergent. 

Let denote the space of rapidly decaying functions Z" C That is, € 

(S(Z") if for any M < oo there exists a constant C^^m such that < C^^m{0~^ 

holds for all ^ e Z". The topology on 5(Z") is given by the seminorms p^, where 

e No and Pk{'^) '■= sup^£2„(0'^ |(/?(^)| . One can show that the continuous linear 
functionals on S{U^) are of the form 

where functions : Z" — > C grow at most polynomially at infinity, i.e. there exist 
constants M < oo and C^^m such that 

\u{0\<C^,m{0'' 

holds for all G Z". Such distributions m : Z" — C form the space 5'(Z"). 
Let : <S(M"^) — >■ 5(M'*) be the Euclidean Fourier transform defined by 

(=^R"/)(0 := / e---« fix) dx, 

where dx = (27r)~"' dx. Mapping : 5(M") — > 5(M") is a bijection, and its inverse 
is given by 

fix) = [ e-« (^Rn/)(0 de. 

JR" 

As it is well-known, this Fourier transform can be uniquely extended to : 
^ S'iR"") by duality. 
Let : C°°(T") ^ 5(Z") be the toroidal Fourier transform defined by 

(2.1) (^T«/)(0 := / e-^^-« /(x) do;, 

where die = (27r)~'" dx. Then ^i n is a bijection and its inverse ^^r^n '■ <S(Z") — > 
C°°(T") is given by 

(^f„i^)(x)= J]e'-« ^(0, 
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SO that f{x) = Egez-e^"^'^ (^T"/)(0- By dualising .^^^^ : 5(Z") ^ C°°(T"), the 
Fourier transform is extended uniquely to mapping : ©'(T"^) 4S'(Z") by the 
formula 

(2.2) {^T"U, (f) := (m, l o ^~„^(/9), 

where m G 'D'(T"), G 5(Z"), and l is defined by {t o '?/')(x) = ■?/'(— x). One can 
easily check that extension fl2.2p when restricted to C°°(T"), is compatible with 
the definition fl2.ip . Here, inclusion C°°(T") C V'{T") is interpreted by {u,ip) = 
m(x) </?(a;) Sx. 

It can be also easily seen that the Fourier analysis on the A^-inflated torus A^T" = 
(R/27rA^Z)" can be obtained from the Fourier analysis of the standard torus T*^ by 

changes of variables. We can identify the dual group A^T" of A^T" with ;^Z". On 
A^T", we can define the Fourier transform ^^T" ■ C°°(A^T") >S(^Z") by 

JOT" 

where 6:y = (27r)~"' dy. Let us denote gN^x) = g{Nx). Then it is easy to see that 
we have the relation 

(2.3) ^NT^g{v)=N"{^T^9N){Nr]), r] e j^Z\ 

3. Difference calculus and discrete Taylor's theorem 

In this section we develop the discrete calculus which will be needed in the sequel. 
In particular, we will formulate and prove a discrete version of the Taylor expansion 
formula on the lattice Z". 

Let a : Z" ^ C. Let ej G N", {ej)j = 1 and {ej)i = if i ^ j. Define the partial 
difference operator A^^. by 

A5,a(0 :=a(e + e,)-a(0, 

and define 

A° = A"^ ■ ■ ■ A?" 

for a G Nq. It can be easily checked that these difference operators commute, i.e. 
that A|A^ = A^A^ for all multi-indices q;,/3 G Nq. Let us summarize properties of 
differences: 

Proposition 3.1. We have 

Let 0, : Z" — > C. Then we have the Leibnitz formula for differences: 

^n^m) = E it) (a?0(o) (Ar"^) (e + /?). 

I3<a ^'^^ 

Also, the "summation by parts" is given by 

(3.1) Yl (A^^)(o = (-1)1-1 E(Ai"<^)(o m, 
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if 


a 


= 0, 




if 


a 


= 1, 




if 


a 


= 2, 



where {A^.(j)){^) = (f){^) — — cj), with the iterative definition for higher order 
differences. In (13. ip we assume that both series are absolutely convergent. 

Proposition 13. II can be proved by induction and we omit the proof. For 6 E 1^^ and 
a e we define = Oi"'^ ■ ■ ■ where Of = 1 and 

(3.2) ef^'^ = ef\e, -k) = e,{e, - 1) . . . {e, - k). 

Then we have A^^^"^ = a^"'^ 5)(°-7)^ analogy to the Euchdean case when d^O" = 
For 6 > 0, let us denote 

(3.3) It := J2 I^" ■■=- J2 ■ 

0<k<b ~b<k<0 

One can think of /| ■ ■ ■ as a discrete version of the one-dimensional integral ■ ■ - d^; 
in this discrete context, the difference takes the role of the differential operator 

d/de 

In the sequel, we adopt the notational conventions 



t6 rki j-ka-i-. 



and so on. 

Lemma 3.2. If 9 E Z and a G No then 

(3.4) 4V£---4'ri = ^^^"^- 

Proof. We observe simple equalities k^^^ = 1, AfcfcW = i and I^.^kk''^^ = 

from which (13. 4p follows by induction. □ 

We note that Lemma 13.21 can be viewed as a discrete trivial version of the funda- 
mental theorem of calculus: /'(O = /(^) ~ /(O) for smooth enough / : M ^ C 
corresponds to /|Ag/(^) = f{6) — /(O) for / : Z ^ C. Lemma 13.21 immediately 
implies its multidimensional version: 

Corollary 3.3. //0 G Z" and a G then 

(3-5) n4a)<SJ--^£:;r"i = ^«'"'. 

where ll]=i h hh ■■■In, where I, := ' ' ' ^fi'^)''^- 

Now we give the multidimensional version of the Taylor expansion formula. 
Theorem 3.4 (Discrete Taylor series on Z"). Let p : Z" — C. Then we can write 

\a\<M 
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with remainder TM^i, 6) satisfying 

(3.6) \/\t^M{U)\<CM |^^"^| max | AfX^ + z/)], 

' ^ ' ' ' \a\=M, ueQ{e) ' ^ ' 

where Q{0) := {z/ G Z" : \uj\ < \9j\ for all j = 1, . . . ,n}. 

Proof. For 7^ a G Nq, let us denote rUa '■= min{j : aj 7^ 0}. For 6 G Z" and 
i G {1, . . . , n}, let us define z/(6', i, k) G Z" by 

z/(^^,^,A;) := (0i,...,0,_i,A;,O,...,O), 

i.e. 

ifl<j<i, 
u{9,i,k)j = Ik, if j = 

I 0, if i < j < n. 
We claim that the remainder can be written in the form 

(3.7) rM{C,9)= Yl ^"(^'^)' 

\a\=M 

where for each a, we have 

n 

(3.8) r.(e, ^) = n ■ ■ ■ C'2r'^ + '^(^' ^(^- 

recall (13. 3p and (13.51) . The proof of (13.81) is by induction. The first remainder term 
ri is of the claimed form, since 

n 
i=l 

where 

here rg. is of the form (13. 8p for a = Ci, m{a) = i and = 1. So suppose that the 
claim (13.81) is true up to order |a| = M. Then 



rM+i(e,^) = TMi^^O)- J2 ^ O^'^^A^iO 

\a\=M 

, I _ n , r V 



\a\=M 



n 



~ 11 ^feOM)^fc(i,2) ■ ■ ■ ^fcO>,) 

|Q|=Mi=i 

b(C + ^^(6*, k{m 

where we used (13.81) and (13. 5p to obtain the last equality. Combining this with the 
equality 

rria 

p(e + 1^(9, ma, k)) - p(0 = if'-^-'^)^ A^'p(e + y{9, 1, £)), 

i=l 
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rAf+i(e,^) -- 



Err tOj tHJA) rHhaj-l) \- ju{i 

11 ^fc(i,l)^fc(i,2) ■ ■ ■ ^k{j,a,) h{i 
\a\=M j=l 



u(e,ma,k{ma ,ama ))i 



1 = 1 



|/3|=M+li=l 

the last step here is just simple tedious book-keeping. Thus the induction proof of 
(13. 8p is complete. Finally, let us prove estimate (13. 6p . By (13. Sp . we obtain 

|A^rM(e,^)| = 



\a\=M 



ll^fc(i,i)^fc(i, 

lQ|=Mi=l 



,2) 



.,(^. Af + m„, A;(m,, a„J)) 
< V ^|^(")| max |A"+X^ + z/)|, 



|a|=A/ 



where in the last step we used (13. 5p . The proof is complete. 



□ 



Remark 3.5. In Theorem 13. 41 we estimated the remainder over the discrete box Q{9), 
but an estimate over a discrete path would have been enough. 

4. PSEUDO-DIFFERENTIAL OPERATORS AND TOROIDAL SYMBOLS 

We note that given a continuous linear operator A : C°°(T") — > C°°(T"'), we can 
consider its toroidal quantization 



where its toroidal symbol a a G C°°(T"' x Z") is uniquely defined by the formula 

where e^{x) := e'"^'^. We note that for pseudo-differential operators this would be just 
another quantization of the same class, see Remark 14. II for a precise statement. 

For C°°(T'"), \eii)®^e C°°(T2") be defined by {ip ® (^)(a;, y) := il){x)^{y). 

If A : C°°(T") ^ ViT^) is a continuous linear operator then one can verify that 

defines the periodic Schwarz distributional kernel Ka G P'(T^"') of operator A G 
£(C°°(T"'), X>'(T")); a common informal notation is 



A(p{x) 



KA{x,y) if{y) Sy. 



T" 
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The convolution kernel /ca € X>'(T^") of A is related to the Schwartz kernel by 
KA_{x,y) = kA{x,x — y), i.e. we have 



A(p{x) = / kAix.x -y) (p{y) &y 
in the sense of distributions. Notice that if kA{x){y) = kA{x,y) then 

^Tn{kAix))iO=(^Aix,0- 

Let mGM, 0<5, p<l. Then the toroidal symbol class S'^^(T" x Z") consists of 
those functions a{x,C,) which are smooth in x for all ^ G Z", and which satisfy 

for every x G T", for every G and for all ^ G Z". The class S^q{T'' x Z") 
will be often denoted by writing simply 5'™(T" x Z"). 

If a G S'^sO^'' X Z'^), we denote by a{X,D) G Op5;;;5(T'^ x Z") the corresponding 
toroidal pseudo-differential operator defined by 

(4.1) a{X,D)f{x)=Y, [ e'^-^-yy^ a{x,0 f{y)^y- 

Naturally, aa[x,D){x,C) = a{x,^). 

To contrast this with Euclidean (Hormander's) symbol classes, we will write b G 
(M" X M") if 6 G C"^(M" X M") and if 

for every x G M", for every a,/? G Ng, and for all ^ G M". If in addition is 27r- 

periodic for every ^ we will write b G S'™^(T" x M"). The corresponding (Euclidean) 
pseudo-differential operator is then given by 

6(X, D)f{x) = [ [ e^(^-^)-« 6(x, fiv) Sy d^ 

The class S'™^(T" x T" x Z") of toroidal compound symbols consists of the functions 
a{x,y,^) which are smooth in x and y for all ^ G Z" and which satisfy 

(4.2) I A^°9f y, I < Caap,m 

for every x,y & T", for every a, /5, 7 G Nq , and for all ^ G Z". Such a function a will 
be also called a toroidal amplitude of order m G M of type {p,S). Formally we may 
also define 

{Op{a)f){x):=J2 I e'(^-^)«a(a;,i/,0/(l/)^2/ 



for / G C°°(T"). 

Remark 4.1. On T", Hormander's usual (p, 5) class of pseudo-differential operators 
Op5'^5(M" X M") of order m G M which are 27r-periodic in x coincides with the class 
Op^m (T" X Z"), i.e. 

OpS'™5(T" X R") = Ops;:;5(T" X Z"), 

see e.g. [16], but in Theorem 15.21 we give a precise relation between actual symbols. 
The relation between the corresponding operators is then given in Theorem 16.41 
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Now we will discuss the calculus of pseudo-differential operators with toroidal sym- 
bols. For this, let us fix the notation first. We define 

K^-^J ^Vn 1 



where Df^ = I and 



Vj Vj 



idyj J idyj \idyj J V%i 



Theorem 4.2 (Compound symbols). Let < 6 < p < 1. For every amplitude 
a E 5'™^(T"' X T" X Z") there exists a unique symbol a E S^^{T"- x Z") satisfying 
Op(a) = 0p((7), where 

(4.4) aix,0 ~ Dl;^aix,y,0\y=.- 

The proof of this theorem is essentially the same as that of the one-dimensional 
case of S']"q(T x Z) analyzed in [32]. In [32] also the following results were proved 
for symbols in S'{"q(T x Z) (see also [25]), and below we give their extensions to 
(p, 6) classes in n dimensions. The possibility for such an extension comes from the 
discrete Taylor expansion and estimates for remainders in Theorem 13. 4[ We omit 
these lengthy proofs and refer to e.g. Theorem 18.21 for a more general version of this 
type of arguments. 

Theorem 4.3 (Composition). Let < 6 < p < 1 and let a a e S'^^{T'' x Z") and 
gb e 5^|(T" X Z"). Then gab e S'^J^"''iT'' x Z") has the following asymptotic 
expansion: 

(4.5) aABix,0 ~ (A^Va(x,0) Di^^aBix,^, 

In particular, the symbol of the commutator [A, B] = AB — BA belongs to the class 
'S'^i^"''"^(T" X Z"). 

Theorem 4.4 (Adjoint). Let < 5 < p < 1. Let ga ^ ^"^(T" x Z"). Then the 
symbol a a* G ^"^(T" x Z") of the adjoint operator A* has the asymptotic expansion 

(4.6) a A* (x, ~ E ^ A°D(")^;(^. 

From the asymptotic expansion for composition of pseudo-differential operators, 
we get an expansion for a parametrix of an elliptic operator: 

Theorem 4.5 (EUipticity) . Let a a G S'™'(T" x Z") be elliptic in the sense that there 
exist constants Co > and Nq E Nq such that 

(4.7) io-^(x,oi>co (er 

for all {x,^) G T" X Z" for which \^\ > Nq; this is equivalent to assuming that there 
exists gb G 5'-™'(T" x Z") such that I — BA, I — AB are smoothing. Let A ~ Yl'jLo 
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where a a, G ^'""^(T" x Z"). Then B ~ J2T=oBk, where Bk E 5-™-'=(T" x Z") ts 
such that (TBg(x, ^) = 1/(T^o(x, ^) for large enough \^\, and recursively 

5. Extension of toroidal symbols 

It is often useful to extend toroidal symbols from T" x Z" to T" x M", ideally getting 
symbols in Hormander's symbol classes. The case of n = 1 and (p, 6) = (1, 0) was 
considered in |32j and [25]. This extension can be done with a suitable convolution 
that respects the symbol inequalities. In the following, 5o,c is the Kronecker delta at 
e Z", i.e. ^0,0 = 1, and 5o,? = if ^ 7^ 0. First we prepare the following useful 
functions G 5(R"): 

Lemma 5.1. There exist functions (pa G iS(R") (for each a E Nq) and a function 
9 G 5(M") such that 

Ve{x) ■= e{x + 2iTk) = 1, 

(^Mn^)|z"(0 = <5o,e and dfi^ur^em = A^>a(0 for all ^ G Z". 

Proof. Let us first consider the 1-dimensional case. Let 6 = di E C°°{M}) such that 

supp(^^i) C (-27r,27r), = ^i(x), 9,{2n - y) + O^iy) = 1 

for X G M and for < y < 27r; these assumptions for 61 are enough for us, and of 
course the choice is not unique. In any case, 9i G (S(M^), so that also ^r6'i G 
If ^ G Z^ then we have 

^^61(0 = / Giix) e-'^-« dx = / (^i(x - 2tt) + 9i{x)) e-'^-« d-x = 5o,^. 

JRI Jo 

If a desired 0q G 5(M^) exists, it must satisfy 

[ e''< d^i^MiO d^= [ e''< = (1 - e^")" [ e'^<UO 

JrI JRi JRI 

due to the bijectivity of ^r : 5(]R"'^) S(M}). Integration by parts leads to the 
formula 

(-ixrei(x) = (l-e'T(^MV«)(x). 

Thus 

f {tE^TOi{x), if < |x| < 2n, 
{^i'(po.)ix) = I ' 1, ifx = 0, 

[ 0, if |x| > 2ti. 

The general n-dimensional case is reduced to the 1-dimensional case, since mapping 
^ = (x ^— s> 6'i(xi)6'i(x2) ■ ■ ■6'i(x„)) G 5(M") has the desired properties. □ 

Theorem 5.2. Let < p < 1 and0<6 <1. Symbol a e S'^gi^'' x Z") is a toroidal 
symbol if and only if there exists an Euclidean symbol a G S^g{T^ x R"') such that 
a = a\jny^^n. Moreover, this extended symbol a is unique modulo S'~°°(T" x M"). 
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The relation between the corresponding pseudo-differential operators will be given 
in Theorem 16.41 

Proof. Let us first prove the "if" part. Let a G ^^^(IR" x M"), and in this part we can 
actually allow any p and 6, for example < p,6 < 1. By the Lagrange Mean Value 
Theorem, if |a| = 1 then 

where i] is on the line from ^ to ^ + a. For a general multi- index a G Nq, we also 
have 

for some rj E Q := + ai] x ■ ■ ■ x [C,n,^n + This can be shown by induction. 

Indeed, let us write a = 5 + 7 for some 6 = ej. Then we can calculate 

for some ( and rj, where we used the induction hypothesis in the first line. Therefore, 
we can estimate 

|A"9fS(x,0| = \d?d^M^,Ok=veQ\ 



Let us now prove the "only if" part. First we show uniqueness. Let a, 6 G ^^^(T" x 
R") and assume that a|T"xZ" = &|t"xZ"- Let c = a — b. Then c G S'^^(T" x M") and 
it satisfies c|t"xZ" =0. If ^ G M" \ Z", choose rj E that is the nearest point (or 
one of the nearest points) to ^. Then we have the first order Taylor expansion 



c(a;,0 = c{x,ri) + 



a: |a| = l 



where 



J2 r„(x,e,e-r/) i^-vT, 

a: |a|=l 

e,^)= / {i-t)d^c{x,^ + t9) dt. 

Jo 



Hence we have |c(a;,.^)| < C (^)'^~''. Continuing the argument inductively for c and 
its derivatives, we obtain the uniqueness modulo S~°°{T^ x M"). 

Let us now show the existence. Let 6 G iS(]R") be as in the Lemma 15. 1[ Define 

a : T" X M" ^ C by 

It is easy to see that a = a|T"xZ"- Furthermore, we have 



|9?9fa(x,0| 
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< r' 



|m-p|a|+<5|/3| 



where we used the summation by parts formula (13. ip . In the last two lines we also 
used that (pa ^ 5(]R") and also a simple fact that for p > we have + r])^ < (O^iv)^ 
and + vr^ivr^ < (0"^, for all C,ry G W. Thus a e 5™5(T" x R"). 



□ 



6. Periodization of pseudo-differential operators 



In this section we describe the relation between operators with Euclidean and 
toroidal quantizations and between operators corresponding to symbols a{x,$,) and 
a = a\jny^^n, given by the operator of periodization of functions. 

First we state a property of a pseudo-differential operator a{X, D) to map the space 
<S(M") into itself, which will be of importance. The following class will be sufficient 
for our purposes, and the proof is straightforward. 

Proposition 6.1. Let a = a(x,^) G C°°(]R" x M*^) and assume that there exist e > 
and G M such that for every a, (3 there are constants Cap and M{a,(3) such that 
the estimate 

holds for all x,^ E M". Then pseudo-differential operator a{X, D) with symbol a{x,^) 
is continuous from S{W^) to 



Before analyzing the relation between operators with Euclidean and toroidal quan- 
tizations, we will describe the periodization operator that will be of importance for 
such analysis. 



Theorem 6.2 (Periodization). The periodization Vf : M" 
S{W') IS defined by 



C of a function f G 



(6.1) 

Then V : S{ 
have 

(6.2) 

and 

(6.3) 



Vf{x) := J2f{x + 27rk). 

C°°(T") is surjective and ||'P/||li(T") < 

P/(X)=^T-/((^M"/)|Z")(X) 
J2f{x + 2lTk)=Y,^''-^ (^R"/)(0- 



Li(iRn). Moreover, we 



In particular, for x = this is the Poisson summation formula. 
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We note that by Theorem 16.21 we may extend the periodization to and also 

this extension is surjective. This is actually rather trivial compared to the smooth 
case of Theorem 16.21 because we can find a preimage / G L^(R") of (7 G L^(T") under 
the periodization mapping V by for example setting / = 5'|[o.27r]" and / = otherwise. 

Proof. Let us show the surjectivity of V : S{W) C~(T"). Let 9 G be a 

function defined in Lemma [5. 1[ Then for any g G C°°(T") it holds that 

V{ge){x) = '^g{x + 2Txk) e{x + 2^k) = g{x) ^ e{x + 2r^k) = g{x), 

where g6 is the product of 6 with 27rZ"-periodic function g on M". We omit the 
straightforward proofs for the other claims. □ 

We note that in the case of the A^-infiated torus A^T" we can use the periodization 
operator Vn instead of V, where Vn '■ 5(M") — > C°°{NT"') can be defined by 

(6.4) (W)(a:) = ^^i.(^R"/|^Z")(a;), xeNT\ 

Let us now establish some basic properties of pseudo-differential operators with 
respect to periodization. 

Definition 6.3. We will say that a function a : M" x R" ^ C is 27r-periodic (we 
will always mean that it is periodic with respect to the first variable x G M") if the 
function x 1-^ a{x,^) is 27rZ"-periodic for all ^. 

As before, we use tilde to denote the restriction of a G C°°(]R" x M") to M" x Z^. 
If a{x,^) is 27r-periodic, we can view it as a function on T*^ x Z", and we write 

Theorem 6.4. Let a = a{x,^) G C°°(M" x M") be 2t: -periodic with respect to x for 
every ^ G M". Assume that for every G Nq there are constants Caf3 and M{a,f3) 
such that estimate 

d:d^^a{x,0 

holds for all x,^ EW^ . Let a = a\jn . Then 

(6.5) Voa{X,D)f = a{X,D)oVf 
for all f G 5(M'^). 

Note that it is not important in this theorem that a is in a symbol class. 

Remark 6.5. Note that by Proposition 16.11 both sides of (16.51) are well-defined func- 
tions in C°°(T"). Moreover, equality (16. 5p can be justified for / in larger classes of 
functions. For example, (16. 5p remains true if / G Cq (M") is a C'' compactly supported 
function for k sufficiently large. 

Proof of Theorem \6.4\ Let / G 5(R"). Then we have 

ViaiX,D)f){x) = 5^ (a(X,D)/)(x + 27rA;) 

fceZ" 

= E / e'^^-+2-'=)«a(x + 27rfc,0 (^M"/)(Ode 



M{a,f3) 
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= ^e^-«a(x,O^T"(P/)(0 
= H(X,D)(P/)(a;). 

As usual, these calculations can be justified in the sense of distributions. □ 

Let us now formulate a useful corollary of Theorem 16.41 that will be of importance 
later, in particular in the proof of Theorem 18.41 If in Theorem 16.41 we take function 
/ such that / = (7|[o,27r]" for some g G C°°(T"), and / = otherwise, it follows 
immediately that g = Vf. Adjusting this argument by shifting the cube [0, 27r]" if 
necessary, we obtain 

Corollary 6.6. Let a = a(x, ^) be as in Theorem "KJ, let g G C°°(T"), and let V be 
an open cube in M" with side length equal to 2tt. Assume that the support of g\y is 
contained in V . Then we have the equality 

a{X,D)g = Voa{X,D){g\v), 

where g\v '■ C is defined as the restriction of g to V and equal to zero outside 

ofV. 

Since we do not always have periodic symbols on M" it may be convenient to 
periodize them. If a{X,D) is a pseudo-differential operator with symbol a(x,^), by 
{Va){X, D) we will denote the pseudo-differential operator with symbol 

(Pa)(a;,0 = a(x + 27rA;, 0- 

This procedure makes sense if, for example, a is in L^(]R"') with respect to x. 

In the following proposition we will assume that supports of symbols and functions 
are contained in [— 7r,7r]". We note that this is not restrictive if these functions are 
already compactly supported. Indeed, if supports of f contained in 

some compact set independent of ^, we can find some G N such that they are 
contained in [—N7t,Ntt]"', and then use the analysis on the A^-inflated torus, with 
periodization operator Vn instead of V, defined in 



Proposition 6.7. Let a = a{x,^) satisfy supp a C [— 7r,7r]" x M" and be such that 
for every G Nq there are constants Cap and M{a,f3) such that the estimate 

holds for all a;, G M". Then we have 

a{X,D)f = {Va){X,D)f + Rf, 



QUANTIZATION OF PSEUDO-DIFFERENTIAL OPERATORS ON THE TORUS 



17 



for all f G C°°(M") with supp / C [— 7r,7r]". Operator R extends to a smoothing 
pseudo- differential operator R : V(R"') — > iS(M"). 

Proof. By the definition we can write 

and let us define i?/ := a(X, D)f— (Va){X, D)f. The assumption on the support of 
a implies that for every x there is only one A; G Z" for which a(x + 27rA;, ^) ^ 0, so the 
sum consists of only one term. It follows that Rf{x) = for x G [— 7r,7r]", because 
for such X the non-zero term corresponds to k = 0. Let now x G M" \ [— vr, vr]". Since 

Rfix) = - I I e^(^-^>«a(x + 27rfc,0/(y)%de 

is just a single term and \x — y\ > Q on supp /, we can integrate by parts with 
respect to ^ any number of times. This implies that R G \['~°° and that Rf decays at 
infinity faster than any power. The proof is complete since the same argument can 
be applied to the derivatives of Rf . □ 

Proposition 16.71 allows us to extend formula of Theorem 16.41 to compact pertur- 
bations of periodic symbols. We will use it when a(X, D) is a sum of a constant 
coefficient operator and an operator with compactly (in x) supported symbol. 

Corollary 6.8. Let a(X, D) he an operator with symbol 

a(x,0 = ai{x,i) + ao(a;,0, 

where ai is as in Theorem \6.4\ ai is 27^ -periodic in x for every ^ G M", and is as 
in Proposition \ 6. % supported in [— tTjTt]" x M". Define 

:= ai(x,0 a; G T",^ G Z". 

Then we have 

(6.6) r oa{X,D)f = b{X,D)oVf + V oRf 

for all f G S{W), and operator R extends to R : V'{W) S{W), so that V o R : 
V'{W) C°°(T"). Moreover, %fai,ao G ^^^(M" x M"), thenbe S'^^iT'' x Z"). 

Remark 6.9. Recalling Remark 16.51 (16.61) can be justified for larger function classes, 
e.g. for / G Cg(]R"') for k sufficiently large (which will be of use in Section [7]). 

Proof of Corollary \6.8{ . By Proposition 16.71 we can write 

a{X, D) = ai(X, D) + {Vao){X, D) + i?, 

with R : V'{W) S{W). Let us define 6(x,0 = ai(x,0 + (7'ao)(a;,0, so that 
b = 6|t"xZ"- The symbol b is 27r-periodic, hence for the operator b{X, D) = ai(X, D) + 
{VaQ){X, D) by Theorem 16.41 we have 

V o b{X, D) = b{X, D)or = ai{X,D)oV + V^o{X, D) o V. 
Since R : V'iW) 5(R"), we also have V o R : V'iW) C~(T"). Finally, 
if ai,ao G ^^^(M'^ x W), then b G ^^^(T" x Z") by Theorem O The proof is 
complete. □ 
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7. Toroidal wave front sets 

Here we shall briefly study microlocal analysis not on the cotangent bundle of the 
torus but on T"^ x Z'^, which is better suited for Fourier series representations. Let 
us define mappings 

(7.1) TTMn : \ {0} ^ §"-\ 7rMn(0 = ^/l^l, 

(7.2) TTTnxRn : X (M" \ {0}) ^ T" x §"-i, njn^^^x,^ = ix,^/\^\). 

We say that K cW \ {0} is a cone in if ^ G K and A > imply G K. 
We say that F C Z" \ {0} is a discrete cone if F = fl for some cone K in M"; 
moreover, if here K is open then F is called and open discrete cone. We omit the 
straightforward proof of the following result: 

Proposition 7.1. F C Z'^yjO} is a discrete cone if and only ifT = Z"'n7r]gi(7rKn(F)). 

Let u e P'(T"). The toroidal wave front set WFt(m) c T" x (Z" \ {0}) is defined 
as follows: we say that (a;o,^o) G x (Z" \ {0}) does not belong to WFj'(m) if and 
only if there exist x ^ C°°(T"') and an open discrete cone F C Z" \ {0} such that 
x{xo) 7^ 0, ^0 e F and 

VA^ > < oo G F : \^t^{xu){0\ < Cn{0'^; 

in such a case we say that ^t"(X'^) decays rapidly in F. 

We say that a pseudo-differential operator A G ^'"(T" x Z") = Op5'"(T" x Z") is 
elliptic at the point (xq, ^o) G T" x (Z" \ {0}) if its toroidal symbol : T" x Z" C 
satisfies 

Wa{xo,0\>C iO"' 

for some constant C > as |,^| — > oo, where ^ G F and F C Z" \ {0} is an open 
discrete cone containing ^o- Should ^ t— (Ta{xo,^) be rapidly decaying in an open 
discrete cone containing then A is said to be smoothing at (xc^o)- The toroidal 
characteristic set of A G \Ef™'(T" x Z") is 

charT(A) := {(xo,^o) G T" x (Z" \ {0}) : A is not elliptic at (xo,^o)}, 

and the toroidal wave front set of A is 

WFr(A) := {(xo,^o) e T" x (Z" \ {0}) : A is not smoothing at (xo,^o)}- 

Now, 

WFt(A) U charT(A) = x (Z" \ {0}) 

because if (x, ^) ^ charj'(A), it means that A is elliptic at (x, ^), and hence not 
smoothing. It can be shown that WFt{A) = if and only if A is smoothing, i.e. 
maps D'(T") to C°°(T"') (equivalently, the Schwartz kernel is C°°-smooth). 

Proposition 7.2. Let A, B e Op^™(T" x Z"). Then 

WFt{AB) C WFt(A) n WFt(5). 

Proof. By Theorem 14.31 applied to pseudo-differential operators A and B we notice 
that the toroidal symbol of AB G Op5'^™(T" x Z") has an asymptotic expansion 

aAB{x,0 ~ (A-a^(a;,0) ^i"^^B(x,0 

a>0 



QUANTIZATION OF PSEUDO-DIFFERENTIAL OPERATORS ON THE TORUS 



19 



where in the latter expansion we have used smooth extensions of toroidal symbols. 
This expansion tells that AB is smoothing at (xcCo) if A or i? is smoothing at 



The notion of the toroidal wave front set is compatible with the action of pseudo- 
differential operators: 

Proposition 7.3. Let u e P'(T") and A e Op5^5(T" x Z"), where < p < 1, 
0<6 <1. Then 

WFt{Au) C WFt(m). 
Especially, if E C°°(T") does not vanish, then WFt{(Pu) = WFt{u). 

We omit the proof which is straightforward. 

We will not pursue the complete analysis of toroidal wave front sets here because 
most of their properties can be obtained from the known properties of the usual wave 
front sets and the following relation, where WF{u) stands for the usual Hormander's 
wave front set of a distribution u. 

Theorem 7.4. Let u G Then 



Proof. Without loss of generality, let u G C''(T"-) for some large k, and let x ^ 
C^(R"') such that supp(x) C (0,271)". If ,^^n(xu) decays rapidly in an open cone 
K C M."- then J^fn{V{xu)) = ^ir"(x'w)|z" decays rapidly in the open discrete cone 
Z" n K. Hence WFt(m) C (T" x Z") n WF(m). 

Next, we need to show that (T" x Z") \ WFt{u) C (T" x Z") \ WF(n). Let 
(xo,^o) e (J"" X Z") \ WFt(«) (where 7^ 0). We must show that (a;o,^o) ^ WF(u). 
There exist x ^ C°°(T'"') (we may assume that supp(x) C (0,27r)"' as above) and an 
open cone K C M."- such that x(^o) 0, E n K and that ^T"('P(x^)) decays 
rapidly in Z" fl K. 

Let Ki C M" be an open cone such that ^0 E Ki G K and that the closure 
KIgKU {0}. Take any function w e C°°(§"-i) such that 



Let a G C°°(M"' x M") be independent of x and such that a{x, ^) = u'(^/|^|) whenever 
1^1 > 1. Then a G S^{W x M"). Let a = a|T«xZ", so that a G S°(T" x Z") by 
Theorem 15.21 

By Corollary [631 we have V{x a{X,D)f) = V{x) a{X, D)(Vf) + V{Rf) for all 
Schwartz test functions /, for a smoothing operator R : V{W^) iS(]R"). By Remark 
[Owe also have V{x a{^,D){xu)) = V{x) a{X, D)(V{xu)) + V{R{xu)), where the 
right hand side belongs to C°°(T"), since its Fourier coefficients decay rapidly on the 
whole H". Therefore also V{x ci{X,D){xu)) belongs to C°°(T"). 

Thus X a{X,D){xu) e C^iM""). Let ^ e Ki such that |^| > 1. Then we have 
^M"(a(X,Z})(xn))(0 = w{^/\^\) ^Rn(xn)(0 = ^R"(x«)(0- Thus ^uAxu) decays 
rapidly on Ki. Therefore {xo,^o) does not belong to WF{u). □ 



□ 



WFt(m) = (T" X Z") n WF(m). 
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8. Calculus of Fourier series operators 

In this section we consider analogues of Fourier integral operators on the torus 
T". We will call such operators Fourier series operators and study their composition 
formulae with pseudo-differential operators on the torus. Fourier series operators will 
be defined by the formula 

(8.1) Tu{x):=J2 f e'(<^(^'«)-^«) a(x, My) % 

where a G (7°°(T" x T" x Z") is a toroidal amplitude and is a real- valued phase 
function. We note that if u G C°°(T"), for function Tu to be well-defined on the 
torus we need that the integral (18.11) is 27r-periodic in x. Therefore, by identifying 
functions on T*^ with 27r-periodic functions on R", we will require that the phase 
function : M" x Z" ^ M is such that function x e^''''^^'''^ is 27r-periodic for all 
^ G Z". Then operator T : C°°(T") —>■ !)'(¥") can be justified in the usual way for 
oscillatory integrals. 

Remark 8.1. Assume that function : M" x Z" ^ M is in C'^ with respect to x 
for all ^ G Z". Assume also that function X ^ e''^(^''5) is 27r -periodic for all ^ G Z*^. 
Then functions x i— >■ d^(j){x,^) are 27r-periodic for all ^ G Z*^ and all a G Nq with 
1 < |a| < k. 

Composition formulae of this section can be compared with those obtained by 
Ruzhansky and Sugimoto [IHl [21], but the assumptions on the regularity or bound- 
edness of higher order derivatives of phases and amplitudes are redundant here thanks 
to the fact that ^ G Z" takes only discrete values. 

We recall the notation for the toroidal version of Taylor polynomials and the cor- 
responding derivatives introduced in (13. 2p and (14. 3p . which will be used in the for- 
mulation of the following theorems. We also define 

(8.2) {-DyY''^ = {-Dy,Y'^^^...{-DyJ'^-\ 
where —Df^ = I and 

<-"''''^"' ^ {-i; - -w,{-wr')- ("1; " ■ 

We also recall definition (14.11) of pseudo-differential operators with toroidal symbols. 

Theorem 8.2 (Composition ToP{X, D)). Let : x Z" ^ M be such that function 
X ^ e'<^(^'«) zs 2TT-penodzc for all ^ G Z". Let T : C°°(T") ©'(¥") be defined by 

(8.3) Tu{x):=y^f e^(<^(^'«)-^-«)a(x,i/,Ow(y)%, 
where the toroidal amplitude a G (7°°(T" x T" x Z") satisfies 

for every x,?/ G T", ^ G Z'^ and a,(3 e N^. Let p G 5^(T" x Z") be a toroidal 
symbol and P = p{X, D) the corresponding pseudo-differential operator. Then the 
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composition TP has the form 

TPu{x) = J2[ e'^'^^"'^)-"-^) c(x,2,0 u{z) dz, 

where 

c{x,z,0 = E / e'(^-^)-(^-«) a{x,y,0 p{y,v) Sy 

satisfies 

\d:d^,c{x,z,o\<Co.pmt {^r"-' 

for every x,z E T", ^ G Z" and G Ng. Moreover, we have an asymptotic 
expansion 

c(x,z,0 (-^.)("^ K^,^,0 A^p(^,0] . 

a>0 ^' 

Furthermore, if0<6<p<l,pe ^^^^(T" x Z") and a G S^siT"" x T" x Z"), i/ien 

c G ^;;;5+^(T" X T" X z"). 

Remark 8.3. We note that if T in fl8.3l) is a pseudo-differential operator with phase 
4>{x,^) = x-^ and amphtude a{x, y,^) = ai^x,^) independent of y, then the asymptotic 
expansion formula of two pseudo-differential operators T o P becomes 

cix, ~ V ^ aix, {-D^Y^^A^piz, 0- 

ci>0 

This is another representation for the composition compared to Theorem 14.31 but the 
amplitude of the composition here is a compound symbol. 

Proof. Let us calculate the composition TP: 

TPu{x) = y2 f e'^'^^'"'^^'^'^^ a{x,y,^) Pu{y) dy 

= E / e^(^(^'«)-^-«) a(x, y, E [/^'-'^ 'piy,v)u{z)SzSy 



where 



Denote 6 := r] — C,, so that by the discrete Taylor formula as in Theorem 13. 4^ we 
formally get 

c{x,z,0 ~ E/ e'(^-^)-%(x, E ^ ^^"^ ^M^' ^) 

/Tn Q;l 

Q^l /Tri 

^,>-n -11- 
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a>0 



y=z- 



Now we have to justify the asymptotic expansion. First we take a discrete Taylor 
expansion and using Theorem 13.41 we obtain 

p{y.^ + 0)= V -.e^-^ A-p(y,o + rM(2/,e,^)- 

\uj\<M 

Let Q{9) = {z/ G Z" : < \9j\ for all j = 1, . . . , n} as in Theorem 13.41 Then 

|A^9,^M(y,^,^)| < C {Of m^^\/\-+-dlp{y,i + v)\ 

veQ{e) 

< C (ey^ max {i + z/)^"l°l-^^ 

< max (z/)l^-H-A/| (^)^-H-M 

The corresponding remainder term in the asymptotic expansion of c(x, 2;,^) is 
Rm{x,z,0 = J2 f e'(^-^>%(x,2/,0rM(l/,e,^)ctl/ 

where we integrated by parts exploiting that (/ — Jfy)e^^^~^'''^ = {Oy e^iv-^)-^^ where 
is the Laplacian with respect to y. Thus we get the estimate 



and we take G N so large that this sum over 6 converges. Hence 

|A-9fa2i?M(x,z,o| < c (^)-+^-^. 

This completes the proof of the first part of the theorem. Finally, we assume that 
a G 5*^5 (¥"■ X X Z"^). Then also the terms in the asymptotic expansion and the 
remainder Rm have corresponding decay properties in the ^-differences, leading to 
the amplitude c G S'^g{T'' x x Z""). This completes the proof. □ 

We now formulate the theorem about compositions of operators in the different 
order. 

Theorem 8.4 (Composition P{X, D)oT). Let : R" x Z" ^Rbe such that function 
X ^ e^H^'O 2TT-penodzc for all ^ G Z". Let T : C°°(T") ^ such that 

Tu{x) ■.= J2[ e^^<^(^'«)"^«) a{x,y,0 n{y) Sy, 

where a G C°°(T" x T" x Z") satisfies 

\d:d^ya{x,y,0\<C^Pm (0"^ 



QUANTIZATION OF PSEUDO-DIFFERENTIAL OPERATORS ON THE TORUS 23 

for every x,y & T", ^ G Z" and a, /? G Ng . Assume that for some C > we have 
(8.4) (O<(V.0(x,O)<C(O 
for every x G T", ^ G Z", and that 



1.5) \d^cPix,0\<Ca{0, d:AU{x,0 



< a 



af3 



for every X G T", ^ G Z" and a, (3 e with \f3\ = l. Letp G S^(T"xZ") be a toroidal 
symbol and let P = p{X, D) be the corresponding pseudo-differential operator. Let 
p{x,^) denote an extension ofp{x,^) to a symbol in S'^(T" x M") as given in Theorem 
IKE Then 

PTu{x) =Y [ e'(*("'«)-"-«) c(x,z,0 u{z) dz, 

where 

\d:d^,cix,z,o\<c^^ {^r^' 

for every x, z G T", ^ G Z" and a, /5 G Nq . Moreover, 

(8.6) c{x,z,i) ^ 9,>(x,r/)|,=v.^(.,0 [e'*(^''^'«)a(|/, ^, 0] \y=.. 

where 

^(x, y, := </>(?/, - + - y) ■ V^x, 0- 

Let us make some remarks about quantities appearing in the asymptotic extension 
formula (18.61) . It is geometrically reasonable to evaluate the symbol p(x,^) at the 
real flow generated by the phase function (p of the Fourier series operator T. This is 
the main complication compared with pseudo-differential operators for which we have 
Proposition 17.31 However, although a-priori symbol p is defined only on T" x Z", we 
can still extend it to a symbol p{x,C,) on T" x M" by Theorem 15.21 so that the restric- 
tion d^p{x, rj) |^=Vi,9i(x,o niakes sense. We also note that function \E'(x, y,^) can not be 
in general considered as a function on T" x T" x Z" because it may not be periodic in 
X and y. However, we can still observe that the derivatives dy [e'*^^'^'^^a(?/, 2;, ^)] \y=x 
are periodic in x and z, so all terms in the right hand side of (18. 6p are well-defined 
functions on T" x T" x Z". 

Finally we note that if G ^^^^(M" x M'^), p G ^^^^(T" x Z"), a G S'^giT'' x T" x Z"), 
and < 5 < p < 1, then we also have c G S'™/^(T" x T" x Z"). 

Proof. To simplify the notation, let us drop writing tilde on p, and denote both 
symbols p and p by the same letter p. There should be no confusion since they 
coincide on x Z". Let P = p{X, D). We can write 



PTu{x) = 5"/ e'^^'-y^-'^ p{x,T]) Tu{y) 6.y 

= H I e'^""^)-Xx, 77) I e'^^^y^^'^-'O a{y,z,0 ui^) dz dy 
= e'('^^"'^)-"«) c(x,2,0 m(^) ct^, 
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where 

c(x, z,0=y^ P{x, v) [ e'(<^(^'«)-^(^'«)+(^-^)"') a{y, z, 

Let us fix some x G M", with corresponding equivalence class [x] G T" which we still 
denote by x. Let C be an open cube with side length equal to 27i centred at x. 
Let X = xi^iU) ^ C'°°(T" X T") be such that < x < 1, xi^iV) = 1 ioi \x — y\ < k 
for some sufficiently small k > 0, and such that supp xi^^ ■) H V C V. The last 
condition means that xi^^ ')\v £ C^{V) is supported away from the boundaries of 
the cube V. Let 

cix, z, i) = c^x, z, + c^^x, z, 0, where 

c^(x, z,0=Yl [ e^('^(^'«)-^(^'«)+(^-^)-^) (1 - xix, y)) a{y, z, pix, v) &y, 

c'\x,z,0 = E / e^W^-«)-^(^'«)+(^-^)-'') x{x,y) a{y,z,Op{x,v) % 

1. Estimate on the support of 1 — x- By making a decomposition into cones 
centred at x viewed as a point in R", it follows that we can assume without loss of 
generality that the support of 1 — x is contained in a set where C < \xj — yj\, for 
some I < j < n. In turn, because of the assumption on the support of x{x, ■)\v it 
follows that C < \xj — yj\ < 2tt — C, for some C > 0. Now we are going to apply 
the summation by parts formula (13 .ip to estimate c\x,z,^). First we notice that it 
follows that 

A^^. e'(^-^) '' = e^(^-^)-(''+^i)- e^i^-yyv = QK-^-yyv(^ qK^j-Vj) - l) ^ on supp 

Hence by the summation by parts formula fl3.1l) we get that 

(8.7) Yl ""'^^''^'^ = ( e'^"^-*) - 1)""^' Yl e'^""^)-" A^i'^'^'pix^r]); 



here sum on the right hand side converges absolutely for large enough A^i. On the 
other hand, we can integrate by parts with operator 

1 



' (V,0(y,O)'-i^y0(z/,O' 

where Jify is the Laplace operator with respect to y, and for which we have e"'^'^^''^-* 
Q^Hy^O^ Note that in view of our assumption (18.41) on 0, we have 

l(V,0(2/,O)' - i ^,0(2/, 01 > l(V,0(2/,O)r > C,{0'. 

Therefore, 

c\x,z,0 =Y. I e^('^(^'^)-'^("'«)+"-^)Lf^ { e-'^-'^x 

x(e'(-.-^.)_l)-^^A:^%(a;,r^) {l-x{x,y)) a(y,z,o} % 
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From the properties of amplitudes, we get 

■q^Jn Jvr\{2-K~c>\xj-yj\>c} 

for all N , if we choose large enough N2 and then large enough A^^i. We can easily 
see that similar estimates work for the derivatives of c^, completing the proof on the 
support of 1 — X- 

2. Estimate on the support of x- Extending p G S^{T^ x Z") to a symbol in 
p e 5'^(T" X M") as given in Theorem 15.2^ we will make its usual Taylor expansion at 
r] = Vx(f>{x,^), so that we have 



p{x,r])= 2^ (9^p(x, V^0(x,O) + 

|«|<M 

.8) + 5Z C'a (r? - VM^, OT ra{x, r/ - V,.</)(x, 0), 

\a\=M 

r^ix, i,r]- V,0(a:, 0) = / (1 - t)^'' d^Pi^, tr] + {1 - t) V,0(x, 0) dt. 

Jo 



Then 



|«|<M \a\=AI 



where 



Now using Corollary 16.61 we can calculate 

= d^p{x, V.0(x,O)x 

X /■ /■ e'^^-^)-" [r/ - V.0(x, 0]" e'^^^^'^)-'^^^'-^" y) a{y, z, dy dr/ 
Jr" Jv 

= V.0(x,O) D^"[e^W^'«)-<^(-'5)+(--^)-^^*(-'«)x(x,y) a(y,z,0]|.=.. 

Since x is identically equal to one for y near x, we obtain the asymptotic formula 
f l8.6p . once the remainders Ra are analyzed. 
3. Estimates for the remainder. Let us first write the remainder in the form 



\y=x 



].9) 



X - V^0(x,O)° ^(y) dy, 
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with 

which is a 27r-periodic function of y. Now, we can use Corollary 16.61 to conclude 
that Ra{x,z,S,) in (18.91) is equal to the periodization with respect to x in the form 

=VxSo,{x,z,C), where 

r,(x,e,r7-V,0(x,O)x 



(8.10) 



I" JV 

where we changed the variables hj 6 = rj — \/x4>{x,C,) and where 

^{x, y, ■■= 0(2/, - e) + - 2/) ■ V.0(x, e) 

and 

r,(a;,e,^) = / (1 -t)*'^ 9^>(a;, V.0(x, + t^) dt. 

Since the periodization does not change the orders in z and ^ it is enough to derive 
the required estimates for Sa{x,z,^). 

Let p e C^{W) be such that p{e) = 1 for |^| < e/2, and p{e) = for \e\ > e, for 
some small e > to be chosen later. We decompose 

Saix, z, = Si{x, z, + Si\x, z, 0, where 
(8.11) X r«(x, 9) [ e'*(^*«) y) a{y, z, 0] dy d^^ 

3.1. Estimate for |^| < e(^). For sufficiently small e > 0, for any < t < 1, 
(Vs0(x, ^) + t9) and (^) are equivalent. Indeed, if we use the inequalities 

{z) < 1 + |z| < V2{z), 

we get 

(V.0(x,O + t^) <{C2V2 + e){0 
^(V.0(x, + t9) >1 + I - l^^l > - l^^l > (Ci - e)(0, 



so we will take e < Ci. This equivalence means that for \9\ < function r^ix, ^, 9) 
is dominated by i(^)^~l"l since p G ^^(T" x M"). We will need two auxiliary estimates. 
The first estimate 
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< 



(5<7 



9 \ 



(8.12) 



— — 17— 5| 



follows from the properties of Tq,. Before we state the second estimate, let us analyze 
the structure of dy e'*^^''^''^^. It has at most |a| powers of terms Vy(j){y, — Vx(l){x, $,), 
possibly also multiplied by at most \a\ higher order derivatives dy(f){y,^). The product 
of the terms of the form Vy(p{y,^) — Vx(p{x,^) can be estimated by C{\y — 
The terms containing no difference Vy(f){y, ^) — Vx4'{x,^) are the products of at most 
\a\/2 terms of the type dy(f){y,^), and the product of all such terms can be estimated 

by C(.^)'"'''^. Altogether, we obtain the estimate 



The second auxiliary estimate now is 



(8.13) 



Now we are ready to prove a necessary estimate for S^{x, z, Let 

(1 - iO'^o) 



i + {0'\x~y\ 

where J^e is the Laplace operator with respect to 9. Integrations by parts with Le 
yield 

Qiix-y)-e 



Si{x,z,^) 



" JV 



{l + {0'\x-y 
9 



2\N 



(1 - iO^e 



r^ix,^,9) D'^y [ e^*(^'^'«) x{x,y) a{y,z,i)\ \ dy d^ 

\{x-y)-e 



.jv (i + (erk-y|2) 



N 



|7|<2Ar 



^« [ e^*(^'^'«) x{x, y) a{y, z, 0] (^p (^^^ r,(x, 9)^ | dy d^. 

Using estimates f l8.12p . fl8.13p and the fact that the measure of the support of function 
9 ^ p{9/ (^)) is estimated by (e(,^))", we obtain the estimate 

(l + (0|y-x|)H 



|7|<2Af 

if we choose N large enough, e.g. N > M 



:i + {0\x-y\ 



2\N 



dy 



\a\ 
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Each derivative of ^) with respect to x or ^ gives an extra power of under 

the integral. Integrating by parts, this amounts to taking more y-derivatives, giving 
a higher power of (^). However, this is not a problem if for the estimate for a given 
number of derivatives of the remainder 5*^(0;, 2;, ^), we choose M = \a\ sufficiently 
large. 

3.2. Estimate for \6\ > e(^). Let us define 



u{x, y,^,9) = {x-y)-9 + ^(x, y, ^) = {x - y) ■ (V.0(x, + ^) + 0(2/, - 0- 

From (18.41) and (18.51) we have 

\VyUj\ = I - + V,0 - V,0| < 2^2(1^1 + (0), 
^^■^^^ |V,cu| > 1^1 - I - V.0I > + (1 - Co\x - y\) iO > cm + (0), 

if we choose k < since \x — y\ < k in the support of x in K (recall that we were 
free to choose k > 0). Let us denote 

a^^{x,y,0 ■= e-'*^^'^''^) D^' Qi^i^,y,0_ 

For any u we have an estimate 

(8.15) |9X,(^'2/,OI<C^(0'''', 



because of our assumption (18. 5p that \dy(f){y,C,)\ < Cy(^). For M = |a| > £ we also 
observe that 



(8.16) Kix,^,e)\<C^, \d':aiy,z,0\<Cp{0"'. 



Let us take ^Ly = i\VyUj\ ^ J2j=ii^yj^)^yj- ^e shown by induction that operator 

Ly has the form 

(8.17) = E P^^^= E c.^5.(V,^)'^9>---9j-^, 



\u\<N |At|=2Af 



where = 2N,\6j\ > 1, J2f=i + \iy\ = 2N. It follows from (ES]) and (EUD that 
\Pu,n\ < C{\e\ + iOf^, since for all 6j we have \dl'uj\ < C{\e\ + (0). By the Leibnitz 
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formula we have 



x/^,"[e^*(^'^'«)x(x,y)a(y,z,0] Syde 

71 +72 +73 =a 



71+72 +73=0 

It foUows now from (18.151) and fl8.16p that 

\Si\x,z,i)\<C [ (|^| + (O)^^(O'"'(erd0 
J\e\>e{0/2 

< (^^^^™+l<^l+"~-'^ 

which yields the desired estimate if we take large enough A^. For the derivatives of 
S^\x, z,^), similar to Part 3.1 for S^, we can get extra powers of 9, which can be 
taken care of by choosing large A^. The proof of Theorem 18.41 is now complete. □ 

Remark 8.5. We could also use the following asymptotic expansion for c based on 
the discrete Taylor expansion from Theorem 13.41 

= E;^[^>(^'^)Uv.^(^,oE / 0^-^e'(^-yy'e'''^^'y'^^a{y,z,O dy 

! /Tri 



a>0 



9. Conditions for L^-boundedness 

In this section we will discuss what conditions on the toroidal symbol a guaran- 
tee the L^-boundedness of the corresponding pseudo-differential operator a{X, D) : 
C°°(T"') — i> ©'(T"). In Theorem 19. 2l we will also prove a result on the L^-boundedness 
of Fourier series operators, and in Remark 19.41 we will give its consequence for the 
boundedness in Sobolev spaces. 

Theorem 9.1. Lei a : T" x Z" ^ C he such that 

(9.1) \d^^a{x,i)\<C /or a// (x,0 e T" X Z", 



and all \(3\ <n + l. Then operator a(X, D) extends to a hounded operator on L^(T"''' 
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We note that compared with several well-known theorems on the L^-boundedness 
of pseudo-differential operators (see e.g. Calderon and Vaillancourt [8j, Coifman and 
Meyer [9], Cordes [10]), Theorem 19. II does not require any regularity with respect to 
the ^-variable. In fact, the boundedness of all partial differences of all orders > 1 
with respect to ^ follows automatically from (19.11) . 

Proof. We can write 

a{X,D)f{x) = 5^a(x,0 (^T"/)(Oe^^-^ 

= Yl (^T"a)(r/,0 (^T"/)(0 e'^-(«+'') 



= E E (^T"a)(u;-e,0 (^T"/)(0, 

where (^t"o)(?7,^) is the rj*'^ Fourier coefficient of a(-,^). Here |(^t"o)(?7, ^)| < 
C so that 



\a{X,D)f\\l,„. 



|a(X,D)/(x)p dx 



^ |^Tr"(a(X,D)/)(.;)p 



< 



sup 



5^(^T"a)(^-e,0 (^T"/)(0 



[^T"a)(^-e,OI X 



X 



sup Y l(^T"a)(a;-e,OI ) E l(^Tr"/)(e)f 



< c 



2 

L2(T") 



by the discrete version of Young's inequality. □ 

As an extension of this result to the setting of Fourier series operators, we have 
the following L^-boundedness property. 

Theorem 9.2. Let T : C°°(T") ^ ViT"^) he defined by 

Tu{x) = J2 e''^^'''''^ «(^''^) i'^T"U){k), 

where : x ^ M and a : x Z" ^ C. Assume that function x t— > e"^*^^'^^ is 
2'K-periodic for every ^ G Z*^, and that for all \a\ <2n + l and \(3\ = 1 we have 

(9.2) \dy{x, k)\<C and a^Af0(a;, k) < C for all (x, k) G T" x Z". 
Assume also that 

(9.3) \V^(l){x,k)-V^(l){x,l)\>C\k-l\ for all x eT"", kj e Z". 
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Then T extends to a bounded operator on L^(T"). 

Note that condition fl9.3p is a discrete version of the usual local graph condition for 
Fourier integral operators, necessary for the local L^-boundedness. We also note that 
conditions on the boundedness of higher order differences of phase and amplitude 
would follow automatically from condition fl9.2p . Therefore, this theorem relaxes as- 
sumptions on the behaviour with respect to the dual variable, compared, for example, 
with the corresponding result by Ruzhansky and Sugimoto [20] . 

Proof. Since for m : T*^ C we have ||M||2,2(Tn) = \\.^T"u\\i'2{i"), it is enough to prove 
that operator 

is bounded from ^^(Z") to ^^(T'^). Let us define Siw{x) = e^-^^^''^ a{x, I) w{l), so that 
S = XlieZ" ^i- From the identity 

(w, S*v)f2(^^„-^ = {Sw,v)^2(^jn^ = / ^ e"''^'''''^ a{x,k) w{k) v{x) dx 

we find that the adjoint 5** to 5* is given by 

{S*v){k) = [ e"'*^"''^) ^(MO v{x) dx 

and so we also have 

{S:v){m) = [ e-'<^(^'"^) v{x) dx = 6iUS*v){l). 

It follows that 

SkS*v{x) = e''^'^'"'^^ a{x,k) {S*v){k) 

= [ e^<^(^'^) a{x,k) e-'*^y^'^ v{y) Sy 

Kki{x,y) v{y) dy, 

where Kki{x^ y) = 5ki e''''^*^^''^)"''^*^^''^^! a(x, k) a{y, k). From fl9.2p and compactness of the 
torus it follows that the kernel K^i is bounded and that ||5'jtS';*i;||i2(']i'n) < C^fc/Hf H^^acj-n). 
In particular, we can trivially conclude that for any > we have the estimate 

On the other hand, we have 

{S*iSkw){m) = [ e-'<^(^'')^(^(5fc^i;)(x)dx 

= Sim [ e'['^("'^)-'^("'^)] a(x, k) ^OM) w{k) dx 
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where Kik{m, jj) = dim^kfj. Jjn e^i'i'(^'^)-'t'(^''')^ a{x,k) a{x,l) dx. Now, if 7^ /, inte- 
grating by parts {2n + l)-times with operator j |v "(K^^ff-v "'(/(i^OP ' using the 
periodicity of a and Vx4> (so there are no boundary terms), we get the estimate 

(9-5) |5f^(m,A.)l<Yf|^, 

where we also used that by the discrete Taylor theorem 13.41 the second condition in 
fl^ implies that 

\Vx(f){x,k) -Vx(j){x,l)\ < C\k-l\ for all x G T", kje Z". 

Estimate flQ.Sp implies 

. . (J 

sup^ |Kifc(m,/i)| = \Kik{l,k)\ < |^_^|2n+i > 

and similarly for sup^ J2m^ ^^^^ have 

C 

(9.6) W^l Sk\\i2(^n)^g2(^n) < 



These estimates for norms \\SkSl\\^2(ju)^L^(j,i-^ and H^f S'A;||^2(an)^<^2(an) in 09.41) and 
(19. 6p . respectively, imply the theorem by a modification of Cotlar's lemma given in 
Proposition [131 which we use with H = and g = L2(T"). □ 

The following statement is a modification of the well-known Cotlar's lemma taking 
into account the fact the operators in our application act on functions on different 
Hilbert spaces. The proof follows [281 P- 280] but there is a difference in how we 
estimate operator norms because we cannot immediately replace operator S* by 5**5' 
in the estimates since they act on functions on different spaces. We omit the proof. 

Proposition 9.3 (Variant of Cotlar's lemma). Let Ti^Q he Hilbert spaces. Assume 
that a family of bounded linear operators {Sj : Ti Q}jei^ ^''^d positive constants 
satisfy 

and A = XljeZ'- '^i^) ^ Then the operator S = J2jez-^ "^i satisfies \\S\\^^g < A. 

Remark 9.4 (Sobolev spaces). By using pseudo-differential operator ^s{D) with 
toroidal symbol (1 + e S(q(T'' x Z"), s e R, we can define Sobolev spaces 

//"(¥") on the torus as spaces of all / e P'(T") such that ^s{D)f E L2(T"). Then 
using Theorems 18.21 18.41 and 19. 2[ we obtain the result on the boundedness of Fourier 
series operators on Sobolev spaces. Namely, let T be a Fourier series operator as in 
Theorem 18.41 Then T is bounded on H^(T"') for all s G M. 

10. Applications to hyperbolic equations 

In this section we will briefly discuss how the analysis can be applied to construct 
global parametrices for the hyperbolic equations on the torus and how to embed 
certain problems in MJ^ into the torus. As already mentioned in the introduction, the 
finite propagation speed of singularities in hyperbolic equations allows one to cut-off 
the equation and the Cauchy data for large x for the local analysis of singularities of 
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solutions for bounded times. Then the problem can be embedded into T" in order 
to apply the periodic analysis developed in this paper. One of the advantages of this 
procedure is that since phases and amplitudes now are only evaluated at ^ G Z" one 
can apply this also for problems with low regularity in ^, in particular to problems 
for weakly hyperbolic equations or systems with variable multiplicities. For example, 
if the principal part has constant coefficients then the loss of regularity occurs only 
in ^ so techniques developed in this paper can be applied. 

Let a{X,D) be a pseudo-differential operator with symbol a = a(x, ^) G S'™(M" x 
M") (with some properties to be specified). There is no difference in the subsequent 
argument if a = a{t, x,^) also depends on t. For a function u = u{t, x) of t G M and 
X G M" we write 



aiX,D)uit,x) = / aix,Oi^R-u)it,Oe'^<d^ 

= 11 e<^~yy^ a{x,0u{t,y)6.ydi. 
Jr" Jr^ 

Let u(t, ■) G L^(R") (0 < t < to) be a solution to the hyperbolic problem 

a . 



(10.1) 



i^u{t,x) = a{X, D)u{t,x), 
u{0,x) = f{x), 



where / G L^(R") is compactly supported. 

Assume now that a{X, D) = ai{X, D) + aQ{X, D) where ai(x, ^) is 27r-periodic and 
ao(x,^) is compactly supported in x (assume even that supp ao(-,0 ^ 
simple example is a constant coefficient symbol ai(x,^) = ai{S,)- Let us also assume 
that supp / C [— TT, tt]"'. 

We will now describe a way to periodize problem fllO.ip . According to Proposition 
16.71 we can replace (110.11) by 



(10.2) 



i^M(t, x) = {ai{x, D) + {Vao){X, D))u{t, x) + Ru{t, x), 
u{x,0) = f{x), 



where the symbol ai + Vao is periodic and i? is a smoothing operator. To study 
singularities of (110. ip . it is sufficient to analyze the Cauchy problem 



(10.3) 



i§vit,x) = Mx,D) + (Pao)(X,Z}))t;(t,x), 
v{x,Q) = f{x) 



since by Duhamel's formula we have WF{u — v) = 0. This problem can be transferred 
to the torus. Let w{t, x) = Vv{-, t){x). By Theorem 16.41 it solves the Cauchy problem 



(10.4) 



i^w(t, x) = {ai{x, D) + Vao{X, D))w{t, x), 
w{x, 0) = Vf{x). 



Now, if a G S*^ is of the first order, the calculus constructed in previous sections yields 
the solution in the form 

w{t,x) = TJ{x) = J2 e'"^^*'"''^ b{t,x,k) ^T"(P/)(fc), 
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where ^) and b(t,x,^) satisfy discrete analogues of the eikonal and transport 

equations. Here we note that ^^j^^iV f){k) = {'^w^f){k). We also note that in 
principal the phase (j)(t,x,k) is defined for discrete values of k E Z", so there is no 
issue of regularity, making this representation potentially applicable to low regularity 
problems and weakly hyperbolic equations. 

Let us give a short example. If the symbol ai(x, ^) = cii(0 has constant coefficients 
and belongs to S^, and belongs to S^, we can find that the phase is given by 
(j)(t,x,k) = x-/c + tai(fc). In particular, Va;0(x, fc) = k. Applying a (X, _D) to w(t,x) = 
Ttf{x) and using the composition formula from Theorem 18.41 we obtain 

a{X,D)TJ{x) =Y. I e^«^-^)-^+*'^i('=» c{t,x,k) f{z) dz, 



where 

(10.5) cit.x^k) ^Y.'-T ^^(^'0L=, KK^.x^k). 



a>0 



since function \1/ in Theorem 18.41 vanishes. From this we can find amplitude h from 
the discrete version of the transport equations, details of which we omit here. Fi- 
nally, we note that we can also have an asymptotic expansion for the amplitude h in 
fll0.5l) in terms of discrete differences and corresponding derivatives c?i"^ instead 
of derivatives and , respectively, if we use Remark 18.51 instead of Theorem 18.41 
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